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Abstract 

We re-examine the bounds on supersymmetric particle masses in light of the 
new E821 data on the muon anomalous magnetic moment and the revised theo- 
retical calculations of its hadronic contributions. The current experimental excess 
is either 1.5a or 3.2a depending on whether e + e~ or r-decay data are used in the 
theoretical calculations for the leading order hadronic processes. Neither result 
is compelling evidence for new physics. However, if one interprets the excess as 
coming from supersymmetry, one can obtain upper mass bounds on many of the 
particles of the minimal supersymmetric standard model (MSSM). Within this 
framework we provide a general analysis of the lightest masses as a function of the 
deviation so that future changes in either experimental data or theoretical calcu- 
lations can easily be translated into upper bounds at the desired level of statistical 
significance. In addition, we give specific bounds on sparticle masses in light of 
the latest experimental and theoretical calculations for the MSSM with universal 
slepton masses, with and without universal gaugino masses. 



In February 2001 the Brookhaven E821 experiment [1] reported evidence for a devi- 
ation of the muon magnetic moment, a M , from the Standard Model expectation of about 
2.7cr. Immediately following that announcement appeared a number of papers analyzing 
the reported excess in terms of various forms of new physics, including supersymmetry 
(SUSY) [2, 3, 4, 5]. Shortly thereafter, errors in the theoretical calculation of the mag- 
netic moment within the Standard Model were discovered. In particular, the sign of 
the light-by-light hadronic contribution to was found to be in error [6], shifting the 
theoretical value by roughly 17 x 10~ 10 in the direction of the E821 data. The resulting 
discrepancy between data and theory was then only about 1.6a, leaving little indication 
of new physics. 

Since that initial rise and ebb of interest in a M , there has been progress both theo- 
retically and experimentally. On the theory side, besides the recalculation of the light- 
by-light scattering, there have also appeared new calculations of the other hadronic 
contributions to [7, 8, 9, 10]. On the experimental front, E821 has announced [11] 
an updated measurement of a M using a data set four times larger than that analyzed 
in Ref. [1]. From their combined data sets, they obtain a world average for a M slightly 
higher than their previous measurement, but with significantly smaller errors. These 
various developments have now placed the discrepancy between experiment and theory 
in the range 1.5a to 3.2cr. 

The current deviation of E821's measurement from the Standard Model provides no 
compelling evidence in favor of new physics. However, the attention paid to this process 
over the last year warrants a re-examination of the bounds that can be placed on new 
physics by the current data, in particular, on new SUSY particles. If the small deviation 
in a M is a sign for new physics, then the SUSY explanation is, for many of us, the most 
exciting of the various proposals, since it implies SUSY at a mass scale not far above the 
weak scale. In particular, it implies a light slepton and a light gaugino, though "light" 
can still be as heavy as many hundred GeV. 

This paper plays two roles. First and more importantly, we present general bounds 
on the spectrum of SUSY masses for any deviation of the muon magnetic moment from 
the theoretical Standard Model value. Any changes in the Standard Model calculation 
or to the data itself will not affect this general analysis; when new data or calculations 
appear, revised bounds can simply be read off the plots contained here. Second, we will 
determine specific la upper bounds on SUSY states in light of the most recent E821 
data and the latest theoretical calculations. 

The existence of these bounds will rely on very simple and clearly stated assumptions 
about the SUSY particle spectrum; these assumptions will not include a fine-tunimg 
constraints. 
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1 SUSY and 

The measurement performed by the E821 collaboration is of the muon's anomalous 
magnetic moment 1 , which is to say, the coefficient a M of the non-renormalizable operator 

^a a ^F aP . (1) 



2m M 



Within the Standard Model, a M receives contributions from QED, electroweak and 
hadronic processes, the latter usually separated into contributions from vacuum po- 
larization and light-by- light scattering. While the QED and electroweak contributions 
are well understood, the hadronic calculations are under scrutiny and require experi- 
mental input. For the purposes of this paper we will use the most recent determination 
of the hadronic contributions available. 

The QED and electroweak contributions to are known: 



a QED = 11658470.56(0.29) [13] 
aj w = 15.2 (0.1) [14] 



where we are expressing in units of 10~ 10 . In the hadronic sector, the next-to-leading 
order (NLO) contributions have been known for several years, while the and light-by- 
light (LBL) contributions have been updated in the last year: 

a Had,NLO = _io.o(0.6) [15] 

a Had,LBL = 8 .0(4.0) [6]. 

The calculation of leading order (LO) hadronic contributions requires input from exper- 
iment. Davier et al. [9] have calculated the LO hadronic pieces using e + e~ scattering 
data and alternatively with a combination of e + e~ and r-decay data. They find: 

Had,LO _ / 684.7 (7.0) (no r data) 
^ ~\ 701.9(6.2) (r data). W 

Hagiwara et al. have also calculated the LO hadronic contribution without the r data 
in Ref. [10]: 

a Had,LO = 683.1 (6.2) (nor data) (5) 

which is consistent with the above evaluation. For our estimate of a M without r data 
we will use the average of the two results stated above and the larger error. Adding the 
contributions gives the standard model prediction for a^. 

J 11659167.7(8.1) (no r data) ( 
a »~\ 11 659 185.7 (7.4) (r data). [b) 



1 See, however Ref. [12]. 
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Figure 1: Super symmetric diagrams contributing to a M at one-loop. 



Since the two values are mutually inconsistent, we will not combine them into a single 
prediction of a M . 

The new measurement made by E821 is [11] ajf 21 = 11 659 204 (7) (5) x 1(T 10 yielding 
a world average of: 

a^ x P = 11 659 203 (8) x 10~ 10 (7) 

from which one deduces a discrepancy between the experiment and the Standard Model 
of 

f 35(11) x 10- 10 (no r data) , . 

da » ~ \ 17(11) x 10- 10 (r data) W 

where we have added the theoretical and experimental errors in quadrature. Thus, the 
deviation is either 3.2a or 1.5a, depending on whether or not the r data is used. In the 
analysis that follows, we present results using both values of Sa^ from Eq. (8). 

The SUSY contributions to a M have been known since the early days of SUSY and 
have become more complete with time [16]. In this paper we will follow the notation 
of Ref. [4] which has the advantage of using the standard conventions of Haber and 
Kane [17]; any convention which we do not define here can be found in either of these 
two papers. 

Prior to the revision of the theoretical calculations and the release of the newest E821 
data, there were a number of analyses in the context of SUSY [2, 3, 4, 5]. The present 
authors also presented an analysis of the full minimal supersymmetric standard model 
(MSSM) with and without gaugino unification and neutralino LSP constraints [18]. 



1.1 The diagrams 

In the mass eigenbasis, there are only two one-loop SUSY diagrams which contribute 
to a M , shown in Figure 1. The first has an internal loop of smuons and neutralinos, the 
second a loop of sneutrinos and charginos. But the charginos, neutralinos and even the 
smuons are themselves admixtures of various interaction eigenstates and we can better 
understand the physics involved by working in terms of these interaction diagrams, of 
which there are many more than two. We can easily separate the leading and sub-leading 
diagrams in the interaction eigenbasis by a few simple observations. 

First, the magnetic moment operator is a helicity- flipping interaction. Thus any 
diagram which contributes to a M must involve a helicity flip somewhere along the fermion 
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current. This automatically divides the diagrams into two classes: those with helicity 
flips on the external legs and those with flips on an internal line. For those in the first 
class, the amplitude must scale as m M ; for those in the second, the amplitudes can scale 
instead by msusY, where msusY represents the mass of the internal SUSY fermion (a 
chargino or neutralino). Since msusY ^ "V, it is the latter class that will typically 
dominate the SUSY contribution to a^. Therefore we will restrict further discussion to 
this latter class of diagrams alone. 

Secondly, the interaction of the neutralinos and charginos with the (s)muons and 
sneutrinos occurs either through their higgsino or gaugino components. Thus each ver- 
tex implies a factor of either (the muon Yukawa coupling) or g (the weak and/or 
hypercharge gauge coupling). Given two vertices, the diagrams therefore scale as |/ 2 , 
gy^ or g 2 . In the Standard Model, y^ is smaller than g by roughly 1CT 3 . In the minimal 
SUSY standard model (MSSM) at low tan/?, this ratio is essentially unchanged, but 
because y^ scales as 1/ cos/?, at large tan/? (~ 60) the ratio can be reduced to roughly 
10 _1 . Thus we can safely drop the y 2 contributions from our discussions, but at large 
tan/3 we must preserve the gy^ pieces as well as the g 2 pieces. 2 

The pieces that we will keep are therefore shown in Figures 2. In Fig. 2(a)-(e) are 
shown the five neutralino contributions which scale as g 2 or gy^; in Fig. 2(f) is the only 
chargino contribution, scaling as gy^. The contributions to from the ith neutralino 
and the mth smuon due to each of these component diagrams are found to be: 



5a, 



487T 2 m% 



F 2 (x im ) X ' 



g\N 2 x X m \X m 2 

gig2Ni\N i2 x m iX m2 



^giypNnNaX^ 
[ -±g 2 y»N t2 N l3 X 2 ml 



(BB) 
(WB) 
(HB) 
(BH) 
(WH) 



(9) 



and for the kth chargino and the sneutrino: 



5a„ 



F?(x k )g 2 y IJl U h2 V kl . (WH) 



24vr 2 



ml 



(10) 



TV C 

The matrices N, U and V are defined in the appendix along with the functions F 2 ' . 
A careful comparison to the equations in the appendix will reveal that we have dropped 
a number of complex conjugations in the above expressions; it has been shown previ- 
ously [3, 4] that the SUSY contributions to a M are maximized for real entries in the mass 
matrices and so we will not retain phases in our discussion. 

In many of the previous analyses of the MSSM parameter space, it was found that it 
is the chargino-sneutrino diagram at large tan (3 that can most easily generate values of 
5a^ large enough to explain the observed discrepancy. From this observation, one can 
obtain an upper mass bound on the lightest chargino and the muon sneutrino. However, 
this behavior is not completely generic. For example, Martin and Wells emphasized 



2 Pieces which are dropped from our discussion are still retained in the full numerical calculation. 
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Figure 2: Diagrams contributing to a M in the interaction eigenbasis. 

that the BB neutralino contribution could by itself be large enough to generate the old 
observed excess in a M , and since it has no intrinsic tan/5 dependence, they could explain 
the old excess with tan (3 as low as 3. We can reproduce their result in a simple way 
because the BB contribution has a calculable upper bound at which the smuons mix at 
45°, rrifa < m^ 2 , and Cm^^ with Ni = B. Then 

\ Sa I _ . < gj rrym^ / \ ( 100GeV\ 2 

where we have used the fact that (XnXi 2 ) < \ and < 3 and have included a 7% two- 
loop suppression factor. Though any real model will clearly suppress this contribution 
somewhat, this is still 10 2 times larger than needed experimentally. 

This pure BB scenario is actually an experimental worst-case, particularly for hadron 
colliders. The only sparticles that are required to be light are a single neutralino (which 
is probably 5-like) and a single jl. The neutralino is difficult to produce, and if stable, 
impossible to detect directly. The neutralino could be indirectly observed in the decay of 
the jl as missing energy, but production of a jl at a hadron machine is highly suppressed. 
In the worst of all possible worlds, E821 could be explained by only these two light 
sparticles, with the rest of the SUSY spectrum hiding above a TeV. Further, even the 
"light" sparticles can be too heavy to produce at a 500 GeV linear collider. While this 
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case is in no way generic, it demonstrates that the E821 excess does not provide any 
sort of no-lose theorem (even at la) for the Tevatron, the LHC or the NLC. 

This raises an important experimental question: how many of the MSSM states must 
be "light" in order to explain the E821 data? In the worst-case, it would appear to be 
only two. Even in the more optimistic scenario in which the chargino diagram dominates 
5a^, the answer naively appears to be two: a single chargino and a single sneutrino. In 
this limit, 



g 2yix m^m ei 



A* {CO) - 247r 2 m 2 



?C 



in / m c\ \ /100GeV\ 2 /tan/3 

< 2600 x 10" 10 x ^— - 

VIOOGeV/ V mo J \ 30 



(12) 

where we have bounded |F 2 C | by 10 by assuming rriy < 1 TeV. But this discussion is 
overly simplistic, as we will see. 



1.2 Mass correlations 

There are a total of 9 separate sparticles which can enter the loops in Fig. 1: 1 sneutrino, 
2 smuons, 2 charginos and 4 neutralinos. The mass spectrum of these 9 sparticles is 
determined entirely by 7 parameters in the MSSM: 2 soft slepton masses(m£, itlr), 2 
gaugino masses (Mi, M 2 ), the //-term, a soft trilinear slepton coupling (Ap) and finally 
tan (3. Of these, plays almost no role at all and so we leave it out of our discussions 
(see the appendix). And in some well- motivated SUSY-breaking scenarios, Mi and M 2 
are also correlated. Thus there are either 5 or 6 parameters responsible for setting 9 
sparticle masses. There are clearly non-trivial correlations among the masses which can 
be exploited in setting mass limits on the sparticles. 

First, there are well-known correlations between the chargino and neutralino masses; 
for example, a light charged Q ~ W implies a light neutral Nj ~ W and vice-versa. 

There are also correlations in mixed systems (i.e., the neutralinos, charginos and 
smuons) between the masses of the eigenstates and the size of their mixings. Consider 
the case of the smuons in particular; their mass matrix is given in the appendix. On 
diagonalizing, the left-right smuon mixing angle is given simply by: 

2m u utan/3 . . 

tan20fi~— ff-— f. 13 
M Ml - M\ v ' 

The chargino contribution is maximized for large smuon mixing and large mixing occurs 
when the numerator is of order or greater than the denominator; since the former is 
suppressed by m M , one must compensate by having either a very large //-term in the 
numerator or nearly equal M L and M R in the denominator, both of which have profound 
impacts on the spectrum. 

There is one more correlation/constraint that we feel is natural to impose on the 
MSSM spectrum: slepton mass universality. It is well-known that the most general 
version of the MSSM produces huge flavor-changing neutral currents (FCNCs) unless 
some external order is placed on the MSSM spectrum. By far, the simplest such order is 



6 



for sparticles with the same gauge quantum numbers to be degenerate. This requirement 
is most stringent in the squark sector, but also holds in the slepton sector due to non- 
observation of \i — > e7 and r /ry. In any case, mechanisms which generate degeneracy 
in the squarks usually do so among the sleptons as well. Thus we assume m fL = m^ L = 
ttil and rrif R = m^ R = m R . Then the mass matrix for the stau sector is identical to 
that of the smuons with the replacement m M — > m T in the off-diagonal elements. This 
enhancement of the mixing in the stau sector by m r /m M ~ 17 implies that m^ < m^. 
In particular, if 

M\M\ < m\ fi 2 tan 2 (3 (14) 

then m\ x < and QED will be broken by a stau vev. Given slepton universality, this 
imposes a constraint on the smuon mass matrix: 



M 2 M 2 > (^1 mJ/i 2 tan 2 /3 (15) 



or on the smuon mixing angle: 

M \m T ) M\-M 2 R 

where M LjR are the positive roots of M\ R . While not eliminating the possibility of 
Ojx ~ 45°, this formula shows that a fine-tuning of at least 1 part in 17 is needed to 
obtain 0(1) mixing. We will not apply any kind of fine-tuning criterion to our analysis, 
yet we will find that this slepton mass universality constraint sharply reduces the upper 
bounds on slepton masses which we are able to find in our study of points in MSSM 
parameter space. 

(As an aside, if one assumes slepton mass universality at some SUSY-breaking mes- 
senger scale above the weak scale, Yukawa-induced corrections will break universality by 
driving the stau masses down. This effect would further tighten our bounds on smuon 
masses and mixings.) 

The above discussion has an especially large impact on the worst-case scenario in 
which the BB contributions dominates da^. For generic points in MSSM parameter 
space, one expects that tan 29^ < 1/17 which reduces the size of the BB contribution 
by a factor of 17. As a byproduct, the masses required for explaining the E821 anomaly 
are pushed back toward the experimentally accessible region. 



2 Numerical results 

Now that we have established the basic principle of our analysis, we will carry it out in 
detail. We will concentrate on three basic cases. The first case is the one most often 
considered in the literature: gaugino mass unification. Here one assumes that the weak- 
scale gaugino mass parameters {Mi and M 2 ) are equal at the same scale at which the 
gauge couplings unify. This implies that at the weak scale M x = (5/3)(a 1 /a 2 )M 2 . The 
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second case we consider is identical to the first with the added requirement that the 
lightest SUSY sparticle (LSP) be a neutralino. This requirement is motivated by the 
desire to explain astrophysical dark matter by a stable LSP. Finally we will also consider 
the most general case in which all relevant SUSY parameters are left free independent 
of each other; we will refer to this as the "general MSSM" case. 

The basic methodology is simple: we put down a logarithmic grid on the space of 
MSSM parameters (Mi, M 2 , rriL, ttir and /i) for several choices of tan/3. The grid 
extends from 10 GeV for Mi, M 2 and fx, and from 50 GeV for rriL,R, up to 2 TeV for all 
mass parameters. For the case in which gaugino unification is imposed M 2 is no longer 
a free parameter and our grid contains 10 8 points. For the general MSSM case our grid 
contains 3 x 10 9 points. Only fi > is considered since that maximizes the value of Sa^. 
Finally, for our limits on tan/3 we used an adaptive mesh routine which did a better job 
of maximizing Sa^ over the space of MSSM inputs. By running with grids of varying 
resolutions and offsets we estimate the error on our mass bounds to be less than ±5%. 

2.1 Bounds on the lightest sparticles 

Perhaps the most important information that can be garnered from the E821 data is 
an upper bound on the scale of sparticle masses. In particular, one can place upper 
bounds on the masses of the lightest sparticle (s) as a function of 5a fl . Here we will 
derive bounds on the lightest slepton and chargino, but we will also derive bounds on 
the lightest several sparticles, independent of their identity. 

These bounds on additional light sparticles provide an important lesson. Without 
them there remains the very real possibility that the E821 data is explained by a pair of 
light sparticles and that the remaining SUSY spectrum is out of reach experimentally. 
But our additional bounds will give us some indication not only of where we can find 
SUSY, but also of how much information we might be able to extract about the fun- 
damental parameters of SUSY — the more sparticles we detect and measure, the more 
information we will have for disentangling the soft-breaking sector of the MSSM. 

In Fig. 3 we have shown the upper mass bounds for the lightest four sparticles 
assuming gaugino mass unification. These bounds are not bounds on individual species 
of sparticles (which will come in the next section and always be larger than these bounds) 
but simply bounds on whatever sparticle happens to be lightest. The important points 
to note are: (i) the maximum values of the mass correspond to the largest value of tan /?, 
which is to be expected given dominance of the chargino diagram at large tan/5; (ii) 
the 1(7 limit without (with) the r data requires at least 2 sparticles to lie below roughly 
490 (990) GeV; (in) for the central value of the E821 data, at least 4 sparticles must lie 
below 585 (915) GeV; and (iv) for low values of tan f3 a maximum value of <5a M is reached 
(we will return to this later). 

The same plots could be produced with the additional assumption that the LSP be 
a neutralino, but we will only show the case for the LSP bound, in Fig. 4. In this figure, 
the solid lines correspond to a neutralino LSP, while the dotted lines are for the more 
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Figure 3: Bounds on the masses of the four lightest sparticles as a function of 5a^ for tan f3 = 
3, 10, 30 and 50. These figures assume gaugino unification only. The vertical dotted lines 
represent the la bounds using the r-decay data (5a ^ = 6) or not (8a ^ = 24). 
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Figure 4: Bound on the mass of the LSP as a function of 5a^ for tan (3 = 3, 5, 10, 30 and 
50. The dotted lines assume gaugino unification only while the solid lines require additionally 
that the LSP is a neutralino. The dashed line is the bound in the general MSSM, calculated 
at tan j3 = 50. 

general case discussed above (i.e., they match the lines in the tan/3 = 50 plot of Fig. 3). 
Notice that for Sa^ > 40 x 10" 10 there is little difference between the cases with and 
without a neutralino LSP. At the extreme upper and lower values of tan /3 there is also 
little difference. It is only for the intermediate values of tan /3 that the mass bound shifts 
appreciably; for tan/3 = 10 it comes down by as much as 50 GeV when one imposes a 
neutralino LSP. 

Finally, we consider the most general MSSM case, i.e., without gaugino unification. 
Here the correlations are much less pronounced, but interesting bounds still exist. For 
example the central value of the E821 data still demands at least 3 sparticles below 
525 (770) GeV. In figure 5 we demonstrate this explicitly by plotting the masses of 
the four lightest sparticles for tan/3 = 50 and a wide range of <5a M . (We also plot 
the mass bound on the LSP in Fig. 4 with the label "General MSSM.") We see that 
dropping the gaugino unification requirement has one primary effect: the mass of the 
LSP is significantly increased. This is because the LSP in the unified case is usually 
a Ni ~ B but isn't itself responsible for generating <5a M . In the general case, the LSP 
must participate in Sa^ (otherwise its mass could be arbitrarily large) and so is roughly 
the mass of the second lightest sparticle in the unified case, whether that be a ft or 
C. Otherwise the differences between the more general MSSM and the gaugino unified 
MSSM are small. 

We have summarized all this data on the LSP in Table 1 where we have shown the 
mass bounds using the la limit of the E821 data on the LSP for various tan/? values 
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Figure 5: Bounds on the masses of the four lightest sparticles as a function of 5a „ for tan (3 = 50. 
The dotted lines assume gaugino unification only while the solid lines are for the general MSSM. 



with our various assumptions. The numbers represent the bounds without (with) the 
inclusion of the r decay data. The last line in the table represents an upper bound for 
any model with tan/? < 50: mLSP < 475GeV(lTeV) for the E821 la lower bounds of 
5a M . But perhaps of equal importance are the bounds on the next 2 lightest sparticles 
(the "2LSP" and "3LSP"): Using a M calculated without r-decay data, m 2 LSP < 485 GeV 
and m 3 LSp < 630 GeV; for the value of calculated with the r data, m 2L sp < 1 TeV, 
while m 3 LSP is pushed somewhat above 1 TeV. Note that the bound on the 2LSP is 
essentially identical to that on the LSP. The central value of the E821 anomaly further 
implies a fourth (third) sparticle below 1 TeV even in this most general case. 



Mass 


General 


Gaugino 


+ Dark 


Bound 


MSSM 


Unification 


Matter 


tan (3 = 3 


205 (331) 


140 (230) 


115 (215) 


5 


235 (395) 


170 (280) 


135 (280) 


10 


280 (475) 


215 (350) 


150 (325) 


30 


340 (750) 


305 (580) 


275 (565) 


50 




475 (1000) 




370 (765) 


365 (765) 



Table 1: Upper bounds on the mass (in GeV) of the lightest sparticle for the general MSSM, 
the MSSM with gaugino mass unification, and the MSSM with gaugino mass unification plus 
a neutralino LSP. The entries represent the la bound without (with) the inclusion of the r 
decay data. The boldfaced tan = 50 entries represent upper bounds over all tan (3 < 50. 
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Figure 6: Bounds on the masses of N\ and /2i as a function of 5a^ for various tan (3 with 
gaugino mass unification assumed. Again, the vertical dotted lines indicate la bounds without 
(with) the inclusion of the r decay data. 

2.2 Bounds on the sparticle species 

In the previous subsection, we derived bounds on the lightest sparticles, independent 
of the identity of those sparticles. Another important piece of information that can be 
provided by this analysis is bounds on individual species of sparticles, for example, on 
the charginos or on the smuons. These bounds will of necessity be higher than those 
derived in the previous section, but still provide important information about how and 
where to look for SUSY. In particular, they can help us gauge the likelihood of finding 
SUSY at Run II of the Tevatron or at the LHC. 

There is one complication in obtaining these bounds. At low tan f3 the data is most 
easily explained by the neutralino diagrams and as such there must be at least one 
light smuon and one light neutralino. At larger tan/3 contributions from the chargino 
diagrams dominate, implying a light chargino and sneutrino. However the correlations 
already discussed preserve the bounds on the various species over the whole range of 
tan (3. A bound on implies a bound on , and a bound on mg implies a bound on 
at least one of the m^., and in certain cases (such as gaugino unification), the converses 
may be true as well. 

We have shown in Fig. 6 the mass bounds on fti and Ni under the assumption of gaug- 
ino unification; a plot for Ci/N 2 will appear later in our discussion of Tevatron physics. 
Note that Ni must lie below 500 (950) GeV, even for large tan /3, thanks to the gaugino 
unification condition, while fix can be heavier but must still lie below 915 (1800) GeV at 
la. 

Finally, we can consider the general MSSM without gaugino unification. The results 
are shown schematically in Fig. 7 where the general MSSM bounds (solid lines) are 
shown alongside the gaugino unification bounds (dashed lines). 

We can see from the figure that the bound on the \i\ is essentially identical to that 
in the gaugino unification picture. However the gaugino masses have shifted, and the 
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Figure 7: Bounds on the masses of fix, C\, and N\ as a function of 5a^ for tan/3 = 50, 
with gaugino unification (solid) and in the general MSSM (dashed). The two lines for the \i\ 
essentially overlap. 

reason is no mystery. Once again, the lightest neutralino is no longer a I?-like spectator 
to the magnetic moment, but is a VT-like partner of a participating W^-like chargino. 

The results of these plots for the current discrepancy are summarized as follows. For 
the MSSM with gaugino unification, the lightest neutralino must fall below 500 (950) GeV 
for a la deviation. The lighter smuon must lie below 915 GeV (1.8 TeV). Without the 
T-decay data, the lighter chargino falls below 890 GeV; using the r data, the bound rises 
above 2 TeV. 

For the worst case, the general MSSM, the smuon bound increases by only about 
50 GeV over the gaugino unification limits and the neutralino bounds increase to 890 GeV 
(over ITeV), only slightly higher than the unified case. However, the bounds on the 
lighter chargino jump above 2 TeV. 

2.3 Bounds on tan/5 

The final bound we will investigate using the E821 data is on tan/3. There had been, 
after the appearance of the original E821 data, some discussion in the literature about 
which values of tan/3 were capable of explaining it. In particular, at lower tan/3 there 
is a real suppression in the maximum size of Sa^. In Figure 8 we have shown the 
maximum attainable value of Sa^ as a function of tan/3 in the general MSSM; adding 
the assumption of gaugino unification changes the figure only slightly. 

The limit in Fig. 8 clearly divides into two regions. At 5a ^ > 36 x 10" 10 the char ffino 
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Figure 8: Bounds on tan (3 as a function of <5a M for the general MSSM. The dotted lines 
represent the la limits, without (with) the inclusion of the r decay data. 



contribution dominates and thus 5a^ oc y^, scaling linearly with tan/3. At lower 8a 
however, both neutralino and chargino contributions can be important so it becomes 
possible to generate Sa^ with much smaller values of tan f3 than would be possible from 
the charginos alone. The E821 data does not therefore imply any bound on tan f3 what- 
soever, neither at la nor at the experimental central value. Further reductions in the 
size of the error bars will not change this result, so long as the central value remains at 
or below its current la upper bound. 



2.4 Implications for the Tevatron 

At its simplest level, the measurement of 5a fl , an anomaly in the lepton sector, has little 
impact on the Tevatron, a hadron machine. In particular, the light smuons associated 
with Sa^ cannot be directly produced at the Tevatron, occurring only if heavier non- 
leptonic states are produced which then decay to sleptons. In the calculation of a M , 
the only such sparticles are the neutralinos and charginos. These states can be copi- 
ously produced and in fact form the initial state for the "gold-plated" SUSY trilepton 
signature. 

Of particular interest for the trilepton signature are the masses of the lighter chargino 
(Ci) and 2nd lightest neutralino (iV 2 ). Studies of mSUGRA parameter space indi- 
cate that the sensitivity to the trilepton signature at Run II/III of the Tevatron de- 
pends strongly on the mass of sleptons which can appear in the gaugino decay chains. 
For heavy sleptons, the Tevatron is only sensitive to gaugino masses in the range [20] 
< 130 to 140 GeV for 10 fb" 1 of luminosity and 145 to 155 GeV for 30 fb" 1 , where 
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Figure 9: Mass bounds on C\ and Ni (where ~ ra^) as a function of 5a^ for tan/? = 3, 
5, 10, 30 and 50. The dotted lines represent the lu bounds without (with) the inclusion of the 
r decay data. This figure assumes gaugino unification and a neutralino LSP. 

the quoted ranges take one from low to high tan/3. However, for light sleptons (below 
about 200 GeV) the range is considerably extended, up to gaugino masses around 190 
to 210 GeV. 

It is impossible in the kind of analysis presented here to comment on the expected 
cross-sections for the neutralino-chargino production (for example, there is no informa- 
tion in a M on the masses of the t-channel squarks) but we can examine the mass bounds 
on C\ and AT 2 . In Fig. 9 we have shown just that: the upper bound on the heavier of 
either C\ or N2 as a function of Sa^ for several values of tan j3. 

A few comments are in order on the figure. First, this figure assumes gaugino unifica- 
tion; dropping that assumption can lead to significantly heavier and unequal masses for 
the N2 and C\. Second, we have also assumed a neutralino LSP; this is to be expected 
since the event topology for the trilepton signal assumes a stable, neutralino LSP. Fi- 
nally, on the y-axis is actually plotted , but in every case we examined with gaugino 
unification, the difference in the maximum masses of C\ and N 2 differed by at most a 
few GeV. This is because they are both dominantly wino-like in the unified case and 
thus have masses ~ M 2 . 

From the figure it is clear that one cannot devise a no-lose theorem for the Tevatron 
from current E821 data. However, if tan/3 is small and Sa^ > 35 x 10 -10 (the current 
central value using only e + e~ data) then the gauginos may be within the Tevatron's 
reach. We must also emphasize that these are upper bounds on the sparticle masses and 
in no way represent best fits or preferred values. Thus, for small tan j3, there is reason 
to hope that the Tevatron will be able to probe the gaugino sector in Run II or III; for 
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Figure 10: Minimum number of sparticles directly observable at a y/s = 500 GeV linear collider 
as a function of 5a^. The bars represent tan/3 = 3, 10, 30 and 50. This graph assumes gaugino 
unification but does not include additional sleptons implied by slepton mass universality. 



larger tan (3 there is little information about SUSY at the Tevatron to be garnered from 



2.5 Implications for a Linear Collider 

A consensus has emerged in favor of building a y/s = 500 GeV linear collider, presumably 
a factory for sparticles with masses below 250 GeV. What does the measurement of a M 
tell us with regards to our chances for seeing SUSY at ^Js = 500 GeV? And how many 
sparticles will be actually accessible to such a collider? 

The analysis of the previous section can put a lower bound on the number of ob- 
servable sparticles at a linear collider as a function of Sa^ and tan (3 and we show those 
numbers as a histogram in Fig. 10. In this figure, we have shown the minimum number 
of sparticles with mass below 250 GeV for tan/3 = 3, 10, 30 and 50, assuming gaugino 
unification. In the graph, the thinner bars represent smaller tan (3. As is to be expected, 
the number of light states increases with increasing Sa^ and with decreasing tan (3. How- 
ever note that there are no tan/3 = 3 lines for 5a^ > 40 x 10~ 10 since there is no way to 
explain such large Sa^ values at low tan f3. 

We see from the figure that at low tan/3, there is a ltr "guarantee" that at least 
1 to 4 sparticles will be produced at a 500 GeV linear collider, depending on whether 
the r-decay is used. This counting does not include extra sleptons due to slepton mass 
universality; for example, a light muon sneutrino also implies light tau and electron 
sneutrinos, and likewise for the charged smuon. We see also that for tan (3 > 30 there is 
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no such guarantee that a 500 GeV machine would produce on-shell sparticles; this is not 
to be taken to mean that one should not expect their production, simply that a M cannot 
guarantee it. 

A similar bar graph can also be made for a 1 TeV machine, though we do not show 
it here. However the relevant numbers can be inferred from Fig. 3. 

Note that any particular class of models, such as gravity- or gauge-mediated SUSY- 
breaking, may not saturate the bounds presented here. That is, other constraints may 
rule out all regions of parameter space in which Sa^ exceeds some maximum value. 
However if a model does allow a particular value of Sa^, then pursuant to the conditions 
discussed here, that model must have at least as many light particles as the number 
given above. 

3 Conclusions 

Deviations in the muon anomalous magnetic moment have long been advertised as a 
key hunting ground for indirect signatures of SUSY. However, the current experimental 
excess is too small, given the large theoretical uncertainties, to provide statistically 
significant evidence for SUSY. 

The most recent Standard Model calculations indicate an excess of either 1.5a or 
3.2a. If one were to accept the E821 anomaly as evidence for SUSY (at either of these 
values), then a number of statements can be inferred at the la confidence level: 

• The lightest sparticle must lie below 475 GeV (1 TeV); 

• For models with unified gaugino masses, there must be at least 2 sparticles with 
masses below 585 GeV (1 TeV); 

• There is no lower bound on tan (5. 

Bounds on individual species of sparticles are weaker, usually falling at or above 1 TeV. 

However, these pessimistic bounds are over all tan/5, which means effectively that 
they are the bounds when tan/3 = 50, the maximum value we considered. For low 
tan/3, the bounds on sparticle masses are much smaller. At tan/3 = 3 and with gaugino 
unification and a neutralino LSP, there must be a sparticle below 115 (215) GeV, within 
the range accessible in the very near future. So, while the bounds placed on SUSY by 
are relatively weak when no constraints are placed on the MSSM, constraining the model 
by demanding gaugino unification or low(er) tan f3 can bring down the mass bounds into 
the experimentally interesting region. 
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Appendix 

The supersymmetric contributions to a M are generated by diagrams involving charged 
smuons with neutralinos, and sneutrinos with charginos. The most general form of the 
calculation, including phases, takes the form (we follow Ref. [4]): 



" 16tt 2 ^ 



'"" (I-LI 2 + \nt\ 2 )F^x tm ) + ^Re[nt m n? m ]F?(x im ) 



6a& = 



E 



I 12m? v ™ 

J_?!V_n^i2_ L i^fl|2M r c^^_ 1 _ 2m c k 



Mc L k c§}F 2 C (x k ) 



(17) 



where i — 1,2,3,4, m — 1,2, and k = 1,2 label the neutralino, smuon and chargino 
mass eigenstates respectively, and 

nf m = \f2g-i_NnXm2 + y^N i3 X ml , 

nL = -^(g 2 N i2 + g 1 N il )X* ml -y^N i3 X* m2 , 



r R 
c k 



(18) 



y^i = fi^^v/ V2mw cos (3 is the muon Yukawa coupling, and g\ >2 are the U(l) hypercharge 
and SU(2) gauge couplings. The loop functions depend on the variables Xi m = m^Jrn%, 
Xk = rn\, I m\ and are given by 



F 2 N (x) 
F?{x) 
F 2 °(x) 



(1-x) 4 
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(1 - x) 3 
2 

3 



2(1 -x) 3 



1 - 6x + 3x 2 + 2x 3 - 6x 2 lux 

1 — x 2 + 2x In x 

2 + 3x — 6x 2 + x 3 + 6x In x 
3-4a; + a; 2 + 21nx 



(19) 



For degenerate sparticles (x = 1) the functions are normalized so that F?(l) = F 2 N (1) = 
F^\l) = F 2 C (1) = 1. We can also bound the magnitude of some of these functions; in 
particular \F 2 N (x)\ < 3 while |F 2 c (x)| is unbounded as x — > 0. 
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The neutralino and chargino mass matrices are given by 

/Mi -m z s w cp m z s w sp \ 

M 2 m z c w Cf3 -m z c w s j3 

-m z s w C(3 m z c w Cf3 -fx 

\ m z s w Sf3 -m z c w sp -fx / 



(20) 



and 




where sp = sin/3, = cos/3 and likewise for 9 W . The neutralino mixing matrix JV^ and 
the chargino mixing matrices £4; and 14; satisfy 

U*M d V^ = dmg(m di ,m d2 ). (22) 

The smuon mass matrix is given in the basis as: 

M 2 _( M 2 L m^At - //tan/3) \ 

M A- ^(^-//tan/3) Ml J ^ 

where 

M 2 L = m| + (s^- ^)m|cos2/3 

M^ = m\ - s 2 w m 2 z cos 2/3 (24) 

for soft masses m| and m^; the unitary smuon mixing matrix X mn is defined by 

XM- 2 Xt = diag (ml, 171%). (25) 

We will define a smuon mixing angle 0^ such that X u = cos 9^ and X 12 = sin#^. In our 
numerical calculations we will set = 0. At low tan /3 we have checked that varying 
Ajx makes only slight numerical difference, while at large tan/3 it has no observable effect 
whatsoever. Finally, the muon sneutrino mass is related to the left-handed smuon mass 
parameter by 

m~ = m\ + cos 2/2- (26) 

The leading 2-loop contributions to Sa^ have been calculated [21] and have been 
found to suppress the SUSY contribution by a factor (4a/ir) log(msuSY /m^) ~ 0.07; we 
will include this 7% suppression in all our numerical results. 
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